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Abstract 

It is studied the MIT-BIH Normal Sinus Rhythm Database using a statistical 
technique of analysis, that is based on the Wavelet and Hilbert Transforms. 
With that technique, it was previously found, that there is a collective and 
intrinsic dynamical behavior up to a scale of 64 heartbeats. Now it is shown, 
that using the Biorthogonal wavelet bior3.1 such a behavior reaches the scale 
1024. That result confirms, that the circulatory system is out of equilibrium. 
According to the Statistical Mechanics of Tsallis, and a recent interpretation 
of G. Wilk et al. respect to the non extensive parameter q, the healthy human 
being is characterized by q = 1.7 ± 0.02. 
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1. Introduction. 

The traditional way of measuring the cardiac activity is through the elec- 
trocardiogram (ECG), that records the voltage generated at the heart. The 
study of the ECG is quite complex because it depends on the position of the 
electrodes, the electrical properties of the skin, and the physiological state 
of the human being [lj]. There are several techniques traditionally used in 
order to analyze the heartbeating. Those techniques are so simple as the ones 
measuring the cardiac mean frequency or as complex as those focused on the 
physical interpretation of the ECG signal. The main purpose of those tech- 
niques used in the study of the ECG is to find the physical and mathematical 
properties to be applied for the prognostics of the sickness of the heart. In 
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order to achieve those objectives, it is common the study of the time series 
corresponding to an ECG which elements are the voltages generated at the 
heart. An ECG is obtained sticking electrodes to the skin, amplifying the 
voltage, and then sending the data to a storing place. The right sampling 
frequency satisfies physiological and technical standards [2)], that are applied 
depending of the kind of the electrical signal, the amount of storing memory, 
the level of noise of an ECG signal, and of the scale of study of those data. 
During the years of study of the ECG it has been found, that applying a 
well known topological theorem of Takens-Sauers (T-S) |3) it is feasible to 
analyze the data of an ECG selecting only the maximum -named "R"- of 
each beat. Taking each maximum, it is feasible to rearrange the time series 
into an equivalent one, as a consequence, it is possible to do a faster analysis. 
Through the years it has been found that, that kind of analysis is enoughly 
robust to be used as a predicting and prognosticating tool for several cardiac 
diseases. 

There are several statistical techniques that use the T-S theorem in order 
to discover characteristics that are common to a set of patients 0, 111] ■ One 
of them, based on the Wavelet and Hilbert transforms is of our interest, 
because unveils a hidden behavior of the heart, that is reminiscent of a system 
out of equilibrium. Such a technique of analysis follows a procedure that is 
similar to the one used in the study of some nonlinear mechanical systems 8], 
that is useful to obtain unusual information of the cardiac dynamics 7|. 
The name of that technique is Cumulative Variation of Amplitude Analysis 
(CVAA), and follows the steps: (i) Select a wavelet [lo|: ^ (m) (m is the 
m-th moment, that indicates that x m ■ ty( m \x)dx = 0) in order to 
analyze the data of an ECG: Xi = x(ti),i = 1,2, ...,n, (ii) Select a set of 
temporal scales: s = 1,2, ...,s , where n G K, (iii) Apply a Continuous 
Wavelet Transform to the time series @, E3]: W { = J2j ^ M ~ t j)] ' 

X{, where Wj are the coefficients of the Wavelet Transform, (iv) Apply the 
Hilbert Transform to the time series: hi = Wj + iHj, where H { = H(Wi) 0, 



10] , (v) To obtain the amplitudes Ai = \/Wf + Hf for the elements of the 
new time series, (vi) Normalize the area under the distribution of amplitudes, 
and (vii) Rescale that distribution such as A4 — > Ai ■ P max and -P(A) _ > 
P(v4j)/P max , where P max is the maximum of the normalized distribution. 
A decade ago, the CVAA was successfully applied to a set of time series 
corresponding to healthy people, giving as a result of that study a common 
Gamma distribution: G v (x) = b u+1 x u e~ bx /r(u+ 1), where b = v/x , being x 
the value that maximizes such distribution. The wavelet used in that study 
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is the Gaussian: g n (x) = C n — — e x . The characteristic parameter of the 

(JjJb 

Gamma distribution is v = 1.8 ±0.1 during the day (6 h.), while during the 
night (6 h.) it is v — 1.4 ± 0.1. The fact of a common Gamma distribution 
suggests, that there is an intrinsic dynamics at the heart characteristic of a 
system out of equilibrium jij]. Even more, that collective behavior is typical 
of system that show a phase transition [12| . It is not still clear, the difference 
between the parameters during the day and during the night hours; according 
to that change is due to the circadian cycle that affects the nervous system. 
In this work, with the aim of exploring the behavior of the data according 
to the wavelet chosen, it is studied by the long term the same set of time 
series corresponding to healthy human beings. In order to study the effect 
of each wavelet on the result of CVAA, that one selected [ltjj correspond 
to three quite different families: An Orthogonal (Daubechies), Biorthogonal, 
and Non Orthogonal: The Gaussian function, which is chosen by the authors 
of CVAA. In this study, it is also identified, the temporal scales for the 
uniqueness of the Gamma function. 

2. Analysis of the data. 

The CVAA is used on the data base Physionet: MIT-BIH Normal Sinus 
Rhythm Database (nsrdb) (l3j |. the same that has been used by the authors 
of this technique. It is selected a set of 18 patients: 5 males of 26 to 45 years 
old, and 13 females of ages between 20 and 50. The time series are obtained 
directly of that data base at Physionet of ECGs, recorded at 128 samples 
per second. The period of each ECG is of the order of 24 h. The interbeat 
intervals RR are obtained using a software provided at Physionet. The length 
of each time series is of the size of 10 5 heartbeats. It is applied the CVAA to 
the data using the next wavelets 0: Daubechies: dbl, db2, db3, Biorthogonal: 
bior3.1, bior3.3, bior3.5, and Gaussian: gausl, gaus2, gaus3, following the 
steps suggested at the Section 1. In this study, the moments of the wavelets 
are not higher because the patients at the hospital are supposed to be almost 
in rest during the day, and that they don't use to follow activities that modify 



1 The notation for wavelets is used according to Ref. [10]. For the Orthogonal and Non 
Orthogonal wavelets it is used : wavabrN, where small letters indicate the abbreviation 
for the wavelet "wavabr" and N corresponds to the N — 1-th moment. For the Biorthog- 
onal wavelet the notation is: biorN r .Nd, where N r .Nd corresponds to the moments N r of 
reconstruction and N<j of decomposition. 
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their physiological state in a specific way: At least that hasn't been indicated 
at the data base Physionet. The range of selected scales for our study is 64 
(« 1 min)— 1024 (« lOmin) heartbeats. The goal is to explore the behavior 
at scales higher than 64 |5j, in order to find the limit of fractality of the 
time series of heartbeats, using CVAA and the set of wavelets and scales. 
The long term study of the heartbeating of healthy people starts applying a 



filter to the time series of RR intervals, according to the algorithm at [14 
which avoids the spurious data. After that step, it is followed the procedure 
of CVAA for each of the wavelets and scales already selected. The complete 
analysis of the cardiac time series is done using a software for the analysis of 
ECGs that has been previously described at Ref. [l8j ]. 

3. Results and discussion. 

In figure 1, it is shown an instance of data collapsing of the cardiac data, 
which has been obtained after applying the CVAA. The analysis corresponds 
to the wavelet bior3.1 and the scale 1024. The parameters of the fit are 
v = 1.46 ± 0.03 and x 2 = 0.26 H. In the tables 1-4, it is presented each 
one of the fitting parameters for a wavelet and scale of study. It is possi- 
ble to see at those tables, that as the scale of analysis increases the fitting 
parameters vary up to 80%. The fractality, which is understood as the 
statistical invariance after rescaling, is only seen using the wavelet bior3.1, 
and for it corresponds an almost constant fitting parameter v = 1.43 ± 0.03 
for the scales 64 to 1024. As it is seen at figure 1 and at the table 1, the 
wavelet bior3.1 provides a good data collapse, and numerical fitting for the 
cardiac data. The analyses corresponding to other wavelets show that the 
fitting parameters change mainly inside those scales. Additionally, for the 
wavelets different than the bior3.1, the data collapse is missed starting at the 
scales: [8, 6, 5] x 2 6 , [13, 7] x 2 6 , [3, 3, 5] x 2 6 corresponding to [dbl,db2,db3], 
[bior3.2,bior3.5] and [gausl,gaus2,gaus3] respectively. It is important to com- 



ment that if the data are not filtered according to Ref. [14]], the distributions 
behave in a quite different way, not necessarily as a Gamma function. In 
other words, the data don't collapse into a single distribution; that result is 
independent of the wavelet selected. 

The results obtained during this work show that, it is not necessary to cut 
the time series in shorter pieces e.g. 6 h. -as it is done at Ref. in order 
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to reduce the non stationarities, and to obtain better fits to the parameters 
of the Gamma function. Periods of 6 h. during the day and 6 h. during the 
night doesn't show meaningful difference in v if it is used the wavelet bior3.1 
for CVAA. On the other hand, for a period of the order of 1 h., it is feasible to 
obtain acceptable data collapse but without using the Hilbert Transform [of. 
That suggests that the difference of v during day and night hours is due to 
the wavelet selected, the size of the time series, and not necessarily to the 
non stationary phenomena at an ECG. 

The integral of the Gaussian wavelet, / = G^ a \x)dx , isn't null as 
corresponds by definition, but this fact apparently doesn't affect a CVAA 
at short scales 0]. According to the results shown here, / ~ is important 
for changing the value of the fitting parameters -mainly v—, for the scales 



shown at tables 1-4. Ref. [15|, |16[ suggests that the physiological state of 
a patient is not supposed to show a change at the distributions after the 
normalization. However, one of the results shown at Ref. j^j seems to con- 
tradict that sentence. Here, the selection of the wavelet used in CVAA is 
fundamental at the time of studying long periods of an ECG. The wavelet 
Biorthogonal, moments 3,5 and the Orthogonal Daubechies, moments 0-2 let 
to obtain a good data collapse, and an acceptable fit to the data for scales 
higher than 64. Those graphical and numerical results are missed slower for 
the Biorthogonal and Orthogonal wavelets, than for the Gaussian one, as the 
scale of analysis increases. The difference of irregularity of the time series of 
cardiac data and that of wavelets is also important. 

The Biorthogonal wavelets are frequently used at the wide field of analysis 



of images and signals [17]. However, it is not common to hear about the 
application of a biorthogonal wavelet for the analysis of ECGs. For the 
analysis of cardiac data, the intrinsic fractality of bior3.1 and bior3.3 (this can 
be tested graphically, and statistically using Detrended Fluctuation Analysis 
(DFA) [15|) is decisive at scales higher than 64 heartbeats. According to 
Ref. [15j, a time series of interbeat intervals of length ~ 10 5 corresponding to 
healthy people shows fractality considering scales of the order of 10 4 in DFA 
and 2 6 in CVAA. The results in this document indicate that the data collapse 
obtained using CVAA holds for scales higher than 1024, of the order of 2 x 10 4 
heartbeats, which is ~ 20% of the data. The fitting parameters at so high 
scales increases slowly as much as 0.3 units. That fractal property of the 
Biorthogonal wavelets used in our study, suggests the answer to questioning 
of why data collapse is almost invariant to such high scales discussed here. At 
the tables 1-4, the fitting parameters corresponding to the wavelet bior3.1, 
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are almost the same up to the highest scale; graphically the data collapse 
is hold on for the range 64-1024. For other wavelets that situation doesn't 
happen: As already commented, the fitting parameters change numerically 
and graphically the data collapsing, and is just hold on functionally, because 
of the data collapsing is missed dramatically for scales higher than 64. 

The long term and long scale study of the ECG of healthy humans ap- 
plying the CVAA: The biorthogonal wavelet bior3.1, suggests that does exist 
an intrinsic behavior of the cardiac dynamics, that is still universaler than 
the already discovered at Ref. |5(. According to this result, the intrinsic 
dynamics of a healthy heart isn't affected by any of the possible non sta- 
tionarities produced during normal activities of a person during a complete 
day, including the perturbation of the circadian rhythm, letting us to un- 
veil a deeper dynamics at the heart, that is commonly masked by a set of 
noises of a different kind. The result described here, generalizes that found 
at Ref. p], but contradicts it in the sense that there is a wavelet for which 
the day and night phases under CVAA are not important for the dynamics of 
the heart. The normalized Gaussian distribution -that characterizes such a 
behavior- is the solution to the Fokker-Planck (FP) equation, that describes 
the probability of finding a particle in a fluid, that moves following a dy- 
namics represented by the Langevin equation. The FP equation represents 
classically a system out of equilibrium. Very recently, it has been found, that 
the solution to the FP is able to be parametrized according to the Statis- 



tics of Tsallis [19] . According to those studies [20] , the Gamma distribution 
is : G q (x) = b^ +1 ^x u (q)e- b ^ x /T[v(q) + 1], where v(q) = (q - l)" 1 being 
1 < q < 2. It is known that when q — y 1, the system is most in equilibrium 
while when q — > 2, the system isn't. Taking into account the results of the 



para metrical study at Ref. [20], substituting v = 1.43 ± 0.02 in the last 
equality, it is found that q = 1.7 ± 0.02. That value of the non extensive 
parameter corresponding to the heartbeating of healthy patients, indicates 
that the heart behaves like a system out of equilibrium -like those repre- 
sented by the FP equation-, and that kind of nonlinear phenomena is able 
to be represented by the Non extensive Statistical Mechanics of Tsallis. 

4. Conclusions. 

Using the CVAA: Wavelet bior3.1 it is found, that the intrinsically col- 
lective behavior of the heartbeating of healthy people can be represented 
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according to the Statistical Mechanics of Tsallis, for which corresponds a 
non extensive parameter q = 1.7 ± 0.02. 

Acknowledgments. 

The author thanks to twice organisms: Promep for academical support, 
and to Conacyt, because of the Project of my authorship: Sep-2003-C02- 
45007. Part of this work was done at Unacar, Mexico. 

References 

[1] H. Riekkinen and P. Rautaharju, Circulation 53 (1) (1976) 40; C. U. 
Vogel et al, Eur. J. Pharmacol. 60 (2004) 461. 

[2] American Heart Association, Circulation 93 (5) (1996) 1043 [Circulation 
electronic pages: |http : / /cir c . aha j ournals . org /cgi / content / full / 93 /5 / 1 
el043]. 

[3] F. Takens, Detecting Strange Attractors in Turbulence, Dynamical Sys- 
tems and Turbulence, Warwick 1980. Proceedings of a symposium held 
at the University of Warwick, edited by A. Dold and B. Eckman, 
Springer, Berlin (1981); T. Sauer et al, Journal of Statistical Physics 
65 (3,4) (1991). 

[4] T. Makikallio, Oulu, University library, Oulu, Finland (1998); S. M. 
Pikkujamsa et al, Circulation 100 (1999) 393; M. C. Teich et al., Heart 
rate variability: Measures and Models, Nonlinear Biomedical Signal Pro- 
cessing, Vol. II, Dynamics Analysis and Modeling, Ch. 6, pp. 159-213, 
edited by M. Akay. IEEE Press, New York (2001). 

[5] P. Ch. Ivanov et al, Nature 38 (3) (1996) 323. 

[6] Y. Kimura et al, Am. J. Physiol. Heart Circ. Physiol. 275 (1998) 1993. 
[7] P. A. Ritto et al, Physica A 349 (2005) 292. 

[8] M. Feldman, Mechanical Systems and Signal Processing 8 (2) (1994) 
119; S. Braun and M. Feldman, Mechanical Systems and Signal Pro- 
cessing 11 (4) (1997) 611. 

[9] J. S. Murguia y E. Campos-Canton, Rev. Mex. Fis. 52 (2) (2006) 155. 



7 



[10] The MathWorks, Inc., The MathWorks- Wavelet Toolbox- Analyze and 
synthetize signals and images using wavelet techniques 1994-2006 
[http: / / www.mathworks.com/products/waveletl . 



[11] J. F. Valdes, Bol. Soc. Mex. Ffe. 16 (2) (2002) 85. 

[12] B. J. West, V. Bhargava y A. L. Goldberger, J. Appl. Physiol. 60 (3) 
(1986) 1089. 

[13] A. L. Goldberger et al, Circulation 101 (23) (2000) e215 [Circulation 



electronic pages: http :/ /circ . aha j ournals . org /c gi / content / full / 1 1 / 23 1 
/e215]. 

[14] K. Ho et al, Circulation 96 (1997) 842. 

[15] S. Havlin et al, Physica A 273 (1999) 46. 

[16] A. A. Aghili et al, Phys. Rev. Lett. 74 (1995) 1254. 

[17] L. Menglong et al, Chemistry Online 65 (2002) 02071 
|http://htxb.icas.ac.cn/2002/c02071.htm] ; A. Kourzi, 13th Euro- 
pean Signal Processing Conference, Antalya, Turkey (2005). 
| http://www.ee.bilkent.edu.tr/~signal/defevent/html/abstract/al469.htm ; 



S. Rein et al, IEEE International Conference on Acoustics, Speech, 
and Signal Processing, Montreal, Canada, IV (2004) 341; Ma et al, J. 
of Zhejiang University SCIENCE A 7 (3) (2006) 361. 

[18] P.A. Ritto, Acalan: Revista de la Universidad Autonoma del Carmen 
40 (2006) 17. 

[19] C. Tsallis, J. Stat. Phys. 52 (1-2) (1988) 479. 

[20] G. Wilk and Z. Wlodarczyk, Phys. Rev. Lett. 84 (13) (2000) 2770; T. 
S. Biro y A. Jakovak (2002); A. N. Petridis, Central European Journal 
of Physics 4 (3) (2006) 363. 



8 



Table 1: After the CVAA, these are the parameters that correspond to the fit of the Gamma function to the distribution of 
amplitudes. The scales shown here are (1 — 4) x 2 6 heartbeats. 





M 


WAVELET 


nx2 6 


Para] 


Daubechies 


Biorthogonal 


Gaussian 




dbl 


db2 


db3 


bior3.1 


bior3.3 


bior3.5 


gausl 


gaus2 


gaus3 




V 


1.18 ±0.04 


1.16 ±0.04 


1.18 ±0.04 


1.41 ±0.03 


1.24 ±0.03 


1.22 ±0.03 


1.08 ±0.06 


0.96 ±0.05 


0.99 ±0.04 


1 


b 


3.04 ±0.07 


3.02 ±0.07 


3.03 ±0.07 


3.24 ±0.05 


3.08 ±0.05 


3.07 ±0.06 


2.95 ±0.12 


2.78 ±0.10 


2.80 ±0.09 




x 2 


0.36 


0.35 


0.37 


0.22 


0.27 


0.32 


0.65 


0.57 


0.53 




V 


1.06 ±0.05 


0.97 ±0.05 


0.99 ±0.05 


1.41 ±0.03 


1.18 ±0.03 


1.10 ±0.04 


1.13 ±0.08 


1.00 ±0.06 


1.01 ±0.06 


2 


b 


2.92 ±0.09 


2.78 ±0.09 


2.78 ±0.10 


3.22 ±0.06 


3.00 ±0.06 


2.91 ±0.08 


2.99 ±0.14 


2.85 ±0.12 


2.84 ±0.11 




x 2 


0.51 


0.57 


0.61 


0.25 


0.33 


0.46 


0.78 


0.74 


0.66 




V 


1.08 ±0.06 


0.94 ±0.05 


0.97 ±0.06 


1.37 ±0.03 


1.07 ±0.03 


0.99 ±0.04 


1.11 ±0.07 


1.08 ±0.08 


1.10 ±0.08 


3 


b 


2.94 ±0.11 


2.77 ±0.10 


2.77 ±0.11 


3.16 ±0.05 


2.86 ±0.06 


2.77 ±0.08 


2.91 ±0.13 


2.90 ±0.15 


2.93 ±0.14 




x 2 


0.62 


0.63 


0.68 


0.25 


0.34 


0.51 


0.76 


0.89 


0.83 




V 


1.11 ±0.07 


0.96 ±0.06 


0.98 ±0.06 


1.40 ±0.03 


1.12 ±0.04 


1.07 ±0.06 


1.15 ±0.07 


1.02 ±0.07 


1.07 ±0.06 


4 


b 


2.98 ±0.13 


2.78 ±0.11 


2.77 ±0.12 


3.19 ±0.05 


2.92 ±0.08 


2.88 ±0.11 


2.93 ±0.13 


2.81 ±0.13 


2.84 ±0.12 




x 2 


0.69 


0.70 


0.77 


0.50 


0.70 


0.66 


0.83 


0.83 


0.80 



Table 2: After the CVAA, these are the parameters that correspond to the fit of the Gamma function to the distribution of 
amplitudes. The scales shown here are (5 — 8) x 2 6 heartbeats. 





M 


WAVELET 


nx2 6 


Para] 


Daubechies 


Biorthogonal 


Gaussian 




dbl 


db2 


db3 


bior3.1 


bior3.3 


bior3.5 


gausl 


gaus2 


gaus3 




V 


1.13 ±0.07 


0.98 ±0.06 


1.06 ±0.07 


1.39 ±0.03 


1.09 ±0.04 


1.05 ±0.06 


1.20 ±0.08 


1.02 ±0.07 


1.07 ±0.06 


5 


b 


2.99 ±0.14 


2.81 ±0.12 


2.86 ±0.13 


3.18 ±0.05 


2.89 ±0.08 


2.83 ±0.12 


2.96 ±0.15 


2.78 ±0.13 


2.86 ±0.11 




x 2 


0.75 


0.77 


0.83 


0.25 


0.51 


0.75 


0.93 


0.87 


0.74 




V 


1.13 ±0.07 


1.10 ±0.09 


1.14 ±0.09 


1.39 ±0.03 


1.10 ±0.04 


1.06 ±0.06 


1.23 ±0.09 


1.08 ±0.07 


1.05 ±0.07 


6 


b 


2.98 ±0.14 


2.93 ±0.16 


2.94 ±0.16 


3.17 ±0.04 


2.90 ±0.08 


2.84 ±0.11 


2.98 ±0.17 


2.85 ±0.13 


2.81 ±0.12 




x 2 


0.79 


0.94 


0.98 


0.22 


0.49 


0.72 


1.06 


0.87 


0.81 




V 


1.10 ±0.08 


1.09 ±0.09 


1.16 ± 0.10 


1.43 ±0.03 


1.12 ±0.05 


1.14 ±0.09 


1.24 ±0.09 


1.13 ±0.08 


1.13 ±0.07 


7 


b 


2.93 ±0.14 


2.90 ±0.16 


2.94 ±0.17 


3.20 ±0.04 


2.91 ±0.08 


2.93 ±0.15 


2.99 ±0.16 


2.90 ±0.14 


2.89 ±0.13 




x 2 


0.82 


1.00 


1.08 


0.24 


0.51 


0.94 


1.01 


0.95 


0.85 




V 


1.09 ±0.07 


1.07 ±0.09 


1.16 ±0.11 


1.43 ±0.03 


1.12 ±0.05 


1.10 ±0.07 


1.33 ±0.10 


1.13 ±0.08 


1.19 ±0.08 


8 


b 


2.90 ±0.14 


2.85 ±0.17 


2.91 ±0.19 


3.20 ±0.04 


2.90 ±0.09 


2.85 ±0.14 


3.07 ±0.17 


2.88 ±0.15 


2.97 ±0.14 




x 2 


0.84 


1.05 


1.19 


0.46 


0.54 


0.91 


1.05 


1.05 


0.90 



Table 3: After the CVAA, these are the parameters that correspond to the fit of the Gamma function to the distribution of 
amplitudes. The scales shown here are (9 — 12) x 2 6 heartbeats. 





GO 
M 


WAVELET 


nx2 6 


Para] 


Daubechies 


Biorthogonal 


Gaussian 




dbl 


db2 


db3 


bior3.1 


bior3.3 


bior3.5 


gausl 


gaus2 


gaus3 




V 


1.11 ±0.07 


1.05 ±0.09 


1.18 ±0.12 


1.44 ±0.03 


1.10 ±0.05 


1.09 ±0.08 


1.44 ±0.11 


1.12 ±0.08 


1.22 ±0.08 


9 


b 


2.90 ±0.14 


2.83 ±0.17 


2.94 ±0.20 


3.20 ±0.04 


2.85 ±0.09 


2.82 ±0.14 


3.17 ±0.17 


2.85 ±0.14 


2.99 ±0.14 




x 2 


0.82 


1.06 


1.22 


0.23 


0.57 


0.99 


1.05 


0.99 


0.95 




V 


1.10 ±0.07 


1.04 ±0.09 


1.16 ±0.12 


1.45 ±0.03 


1.12 ±0.05 


1.10 ±0.08 


1.53 ±0.11 


1.20 ±0.08 


1.25 ±0.08 


10 


b 


2.88 ±0.14 


2.80 ±0.17 


2.92 ±0.20 


3.21 ±0.04 


2.87 ±0.09 


2.83 ±0.14 


3.24 ±0.18 


2.91 ±0.15 


3.00 ±0.15 




x 2 


0.82 


1.13 


1.22 


0.23 


0.62 


0.97 


1.01 


1.07 


0.98 




V 


1.12 ±0.08 


0.98 ±0.09 


1.06 ±0.10 


1.45 ±0.03 


1.12 ±0.05 


1.12 ±0.09 


1.68 ±0.12 


1.23 ±0.08 


1.22 ±0.08 


11 


b 


2.89 ±0.14 


2.73 ±0.16 


2.81 ±0.19 


3.21 ±0.04 


2.89 ±0.09 


2.87 ±0.16 


3.40 ±0.19 


2.95 ±0.14 


2.97 ±0.14 




x 2 


0.85 


1.11 


1.24 


0.23 


0.63 


1.07 


1.03 


1.06 


0.95 




V 


1.13 ±0.08 


0.98 ±0.09 


1.04 ±0.10 


1.43 ±0.03 


1.15 ±0.06 


1.14 ±0.09 


1.77 ±0.13 


1.29 ±0.09 


1.26 ±0.09 


12 


b 


2.90 ±0.14 


2.71 ±0.17 


2.76 ±0.18 


3.19 ±0.04 


2.93 ±0.11 


2.92 ±0.16 


3.48 ±0.20 


3.01 ±0.15 


3.01 ±0.15 




x 2 


0.87 


1.18 


1.27 


0.24 


0.69 


1.04 


1.05 


1.11 


1.03 



Table 4: After the CVAA, these are the parameters that correspond to the fit of the Gamma function to the distribution of 
amplitudes. The scales shown here are (13 — 16) x 2 6 heartbeats. 





GO 
M 


WAVELET 


nx2 6 


Para] 


Daubechies 


Biorthogonal 


Gaussian 




dbl 


db2 


db3 


bior3.1 


bior3.3 


bior3.5 


gausl 


gaus2 


gaus3 




V 


1.15 ±0.09 


1.07 ±0.11 


1.12 ±0.11 


1.47 ±0.03 


1.18 ±0.06 


1.10 ±0.08 


1.81 ±0.13 


1.34 ±0.09 


1.29 ±0.09 


13 


b 


2.92 ±0.16 


2.79 ±0.19 


2.85 ±0.20 


3.23 ±0.05 


2.96 ±0.11 


2.88 ±0.15 


3.51 ±0.19 


3.04 ±0.16 


3.01 ±0.15 




x 2 


0.93 


1.29 


1.31 


0.26 


0.70 


0.98 


1.03 


1.12 


1.01 




V 


1.16 ±0.09 


1.07 ±0.10 


1.17 ±0.11 


1.45 ±0.03 


1.09 ±0.05 


1.04 ±0.09 


1.86 ±0.13 


1.40 ±0.10 


1.28 ±0.09 


14 


b 


2.93 ±0.16 


2.80 ±0.17 


2.88 ±0.20 


3.21 ±0.04 


2.85 ±0.10 


2.76 ±0.16 


3.55 ±0.20 


3.10 ±0.16 


3.00 ±0.15 




x 2 


0.93 


1.19 


1.32 


0.24 


0.66 


1.09 


1.06 


1.16 


1.06 




V 


1.18 ±0.09 


1.13 ±0.11 


1.26 ±0.13 


1.46 ±0.03 


1.08 ±0.06 


1.08 ±0.10 


1.78 ±0.12 


1.46 ±0.11 


1.27 ±0.09 


15 


b 


2.96 ±0.17 


2.84 ±0.19 


2.99 ±0.22 


3.22 ±0.05 


2.81 ±0.10 


2.77 ±0.17 


3.49 ±0.19 


3.15 ±0.17 


2.95 ±0.16 




x 2 


0.97 


1.28 


1.42 


0.25 


0.69 


1.21 


1.02 


1.18 


1.13 




V 


1.17 ±0.09 


1.14 ±0.12 


1.31 ±0.14 


1.46 ±0.03 


1.09 ±0.06 


1.19 ±0.10 


1.81 ±0.13 


1.64 ±0.12 


1.37 ±0.11 


16 


b 


2.94 ±0.17 


2.87 ±0.20 


3.06 ±0.24 


3.23 ±0.05 


2.83 ±0.11 


2.91 ±0.18 


3.52 ±0.20 


3.31 ±0.19 


3.07 ±0.18 




x 2 


1.00 


1.34 


1.49 


0.52 


0.75 


1.18 


1.06 


1.22 


1.21 




xP(max) 

Figure 1: A distribution of the data (small black squares) corresponding to NSRDB. The 
parameters of the fit are v = 1.46 ± 0.03 and x 2 = 0.26. 
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